Population genetics of Ice Age brown bears

The Pleistocene was a dynamic period for Holarctic mammal species, complicated by episodes of glaciation, local extinctions, and intercontinental migration. The genetic consequences of these events are difficult to resolve from the study of present-day populations. To provide a direct view of population genetics in the late Pleistocene, we measured mitochondrial DNA sequence variation in seven permafrost-preserved brown bear (Ursus arctos) specimens, dated from 14,000 to 42,000 years ago. Approximately 36,000 years ago, the Beringian brown bear population had a higher genetic diversity than any extant North American population, but by 15,000 years ago genetic diversity appears similar to the modern day. The older, genetically diverse, Beringian population contained sequences from three clades now restricted to local regions within North America, indicating that current phylogeographic patterns may provide misleading data for evolutionary studies and conservation management. The late Pleistocene phylogeographic data also indicate possible colonization routes to areas south of the Cordilleran ice sheet.

The evolution of RNA viruses: A population genetics view

RNA viruses are excellent experimental models for studying evolution under the theoretical framework of population genetics. For a proper justification of this thesis we have introduced some properties of RNA viruses that are relevant for studying evolution. On the other hand, population genetics is a reductionistic theory of evolution. It does not consider or make simplistic assumptions on the transformation laws within and between genotypic and phenotypic spaces. However, such laws are minimized in the case of RNA viruses because the phenotypic space maps onto the genotypic space in a much more linear way than on higher DNA-based organisms. Under experimental conditions, we have tested the role of deleterious and beneficial mutations in the degree of adaptation of vesicular stomatitis virus (VSV), a nonsegmented virus of negative strand. We also have studied how effective population size, initial genetic variability in populations, and environmental heterogeneity shapes the impact of mutations in the evolution of vesicular stomatitis virus. Finally, in an integrative attempt, we discuss pros and cons of the quasispecies theory compared with classic population genetics models for haploid organisms to explain the evolution of RNA viruses.

Multilocus analysis of extracellular putative virulence proteins made by group A Streptococcus: Population genetics, human serologic response, and gene transcription

Species of pathogenic microbes are composed of an array of evolutionarily distinct chromosomal genotypes characterized by diversity in gene content and sequence (allelic variation). The occurrence of substantial genetic diversity has hindered progress in developing a comprehensive understanding of the molecular basis of virulence and new therapeutics such as vaccines. To provide new information that bears on these issues, 11 genes encoding extracellular proteins in the human bacterial pathogen group A Streptococcus identified by analysis of four genomes were studied. Eight of the 11 genes encode proteins with a LPXTG(L) motif that covalently links Gram-positive virulence factors to the bacterial cell surface. Sequence analysis of the 11 genes in 37 geographically and phylogenetically diverse group A Streptococcus strains cultured from patients with different infection types found that recent horizontal gene transfer has contributed substantially to chromosomal diversity. Regions of the inferred proteins likely to interact with the host were identified by molecular population genetic analysis, and Western immunoblot analysis with sera from infected patients confirmed that they were antigenic. Real-time reverse transcriptase–PCR (TaqMan) assays found that transcription of six of the 11 genes was substantially up-regulated in the stationary phase. In addition, transcription of many genes was influenced by the covR and mga trans-acting gene regulatory loci. Multilocus investigation of putative virulence genes by the integrated approach described herein provides an important strategy to aid microbial pathogenesis research and rapidly identify new targets for therapeutics research.

Polymorphic simple sequence repeat regions in chloroplast genomes: applications to the population genetics of pines.

Simple sequence repeats (SSRs), consisting of tandemly repeated multiple copies of mono-, di-, tri-, or tetranucleotide motifs, are ubiquitous in eukaryotic genomes and are frequently used as genetic markers, taking advantage of their length polymorphism. We have examined the polymorphism of such sequences in the chloroplast genomes of plants, by using a PCR-based assay. GenBank searches identified the presence of several (dA)n.(dT)n mononucleotide stretches in chloroplast genomes. A chloroplast (cp) SSR was identified in three pine species (Pinus contorta, Pinus sylvestris, and Pinus thunbergii) 312 bp upstream of the psbA gene. DNA amplification of this repeated region from 11 pine species identified nine length variants. The polymorphic amplified fragments were isolated and the DNA sequence was determined, confirming that the length polymorphism was caused by variation in the length of the repeated region. In the pines, the chloroplast genome is transmitted through pollen and this PCR assay may be used to monitor gene flow in this genus. Analysis of 305 individuals from seven populations of Pinus leucodermis Ant. revealed the presence of four variants with intrapopulational diversities ranging from 0.000 to 0.629 and an average of 0.320. Restriction fragment length polymorphism analysis of cpDNA on the same populations previously failed to detect any variation. Population subdivision based on cpSSR was higher (Gst = 0.22, where Gst is coefficient of gene differentiation) than that revealed in a previous isozyme study (Gst = 0.05). We anticipate that SSR loci within the chloroplast genome should provide a highly informative assay for the analysis of the genetic structure of plant populations.

Social evolution in a new environment: the case of introduced fire ants.

The inadvertent introduction of the fire ant Solenopsis invicta to the United States from South America provides the opportunity to study recent social evolution by comparing social organization in native and introduced populations. We report that several important elements of social organization in multiple-queen nests differ consistently and dramatically between ants in Argentina and the United States. Colonies in Argentina contain relatively few queens and they are close relatives, whereas colonies in the United States contain high numbers of unrelated queens. A corollary of these differences is that workers in the native populations are significantly related to the new queens that they rear in contrast to the zero relatedness between workers and new queens in the introduced populations. The observed differences in queen number and relatedness signal a shift in the breeding biology of the introduced ants that is predicted on the basis of the high population densities in the new range. An additional difference in social organization that we observed, greater proportions of permanently unmated queens in introduced than in native populations, is predicted from the loss of alleles at the sex-determining locus and consequent skewing of operational sex ratios in the colonizing ants. Thus, significant recent social evolution in fire ants is consistent with theoretical expectations based on the altered ecology and population genetics of the introduced populations.

Concordance of gene genealogies reveals reproductive isolation in the pathogenic fungus Coccidioides immitis

Simple cladogenetic theory suggests that gene genealogies can be used to detect mixis in a population and delineate reproductively isolated groups within sexual taxa. We have taken this approach in a study of Coccidioides immitis, an ascomycete fungus responsible for a recent epidemic of coccidioidomycosis (Valley fever) in California. To test whether this fungus represents a single sexual species throughout its entire geographic range, we have compared genealogies from fragments of five nuclear genes. The five genealogies show multiple incompatibilities indicative of sex, but also share a branch that partitions the isolates into two reproductively isolated taxa, one centered in California and the other outside California. We conclude that coccidioidomycosis can be caused by two distinct noninterbreeding taxa. This result should aid the future study of the disease and illustrates the utility of the genealogical approach in population genetics.

Expression of alfalfa mosaic virus coat protein in tobacco mosaic virus (TMV) deficient in the production of its native coat protein supports long-distance movement of a chimeric TMV

Alfalfa mosaic virus (AlMV) coat protein is involved in systemic infection of host plants, and a specific mutation in this gene prevents the virus from moving into the upper uninoculated leaves. The coat protein also is required for different viral functions during early and late infection. To study the role of the coat protein in long-distance movement of AlMV independent of other vital functions during virus infection, we cloned the gene encoding the coat protein of AlMV into a tobacco mosaic virus (TMV)-based vector Av. This vector is deficient in long-distance movement and is limited to locally inoculated leaves because of the lack of native TMV coat protein. Expression of AlMV coat protein, directed by the subgenomic promoter of TMV coat protein in Av, supported systemic infection with the chimeric virus in Nicotiana benthamiana, Nicotiana tabacum MD609, and Spinacia oleracea. The host range of TMV was extended to include spinach as a permissive host. Here we report the alteration of a host range by incorporating genetic determinants from another virus.

Does recombination improve selection on codon usage? Lessons from nematode and fly complete genomes

Understanding the factors responsible for variations in mutation patterns and selection efficacy along chromosomes is a prerequisite for deciphering genome sequences. Population genetics models predict a positive correlation between the efficacy of selection at a given locus and the local rate of recombination because of Hill–Robertson effects. Codon usage is considered one of the most striking examples that support this prediction at the molecular level. In a wide range of species including Caenorhabditis elegans and Drosophila melanogaster, codon usage is essentially shaped by selection acting for translational efficiency. Codon usage bias correlates positively with recombination rate in Drosophila, apparently supporting the hypothesis that selection on codon usage is improved by recombination. Here we present an exhaustive analysis of codon usage in C. elegans and D. melanogaster complete genomes. We show that in both genomes there is a positive correlation between recombination rate and the frequency of optimal codons. However, we demonstrate that in both species, this effect is due to a mutational bias toward G and C bases in regions of high recombination rate, possibly as a direct consequence of the recombination process. The correlation between codon usage bias and recombination rate in these species appears to be essentially determined by recombination-dependent mutational patterns, rather than selective effects. This result highlights that it is necessary to take into account the mutagenic effect of recombination to understand the evolutionary role and impact of recombination.

Molecular markers reveal that population structure of the human pathogen Candida albicans exhibits both clonality and recombination.

The life history of Candida albicans presents an enigma: this species is thought to be exclusively asexual, yet strains show extensive phenotypic variation. To address the population genetics of C. albicans, we developed a genetic typing method for codominant single-locus markers by screening randomly amplified DNA for single-strand conformation polymorphisms. DNA fragments amplified by arbitrary primers were initially screened for single-strand conformation polymorphisms and later sequenced using locus-specific primers. A total of 12 single base mutations and insertions were detected from six out of eight PCR fragments. Patterns of sequence-level polymorphism observed for individual strains detected considerable heterozygosity at the DNA sequence level, supporting the view that most C. albicans strains are diploid. Population genetic analyses of 52 natural isolates from Duke University Medical Center provide evidence for both clonality and recombination in C. albicans. Evidence for clonality is supported by the presence of several overrepresented genotypes, as well as by deviation of genotypic frequencies from random (Hardy-Weinberg) expectations. However, tests for nonrandom association of alleles across loci reveal less evidence for linkage disequilibrium than expected for strictly clonal populations. Although C. albicans populations are primarily clonal, evidence for recombination suggests that sexual reproduction or some other form of genetic exchange occurs in this species.

Cloning of inversion breakpoints in the Anopheles gambiae complex traces a transposable element at the inversion junction

Anopheles arabiensis, one of the two most potent malaria vectors of the gambiae complex, is characterized by the presence of chromosomal paracentric inversions. Elucidation of the nature and the dynamics of these inversions is of paramount importance for the understanding of the population genetics and evolutionary biology of this mosquito and of the impact on malaria epidemiology. We report here the cloning of the breakpoints of the naturally occurring polymorphic inversion 2Rd′ of A. arabiensis. A cDNA clone that cytologically mapped on the proximal breakpoint was the starting material for the isolation of a cosmid clone that spanned the breakpoint. Analysis of the surrounding sequences demonstrated that adjacent to the distal breakpoint lies a repetitive element that exhibits distinct distribution in different A. arabiensis strains. Sequencing analysis of that area revealed elements characteristic of transposable element terminal repeats. We called this presumed transposable element Odysseus. The presence of Odysseus at the junction of the naturally occuring inversion 2Rd′ suggests that the inversion may be the result of the transposable element’s activity. Characteristics of Odysseus’ terminal region as well as its cytological distribution in different strains may indicate a relatively recent activity of Odysseus.

ALFRED: an allele frequency database for diverse populations and DNA polymorphisms—an update

ALFRED (the ALelle FREquency Database) is designed to store and disseminate frequencies of alleles at human polymorphic sites for multiple populations, primarily for the population genetics and molecular anthropology communities. Currently ALFRED has information on over 180 polymorphic sites for more than 70 populations. Since our initial release of the database we have focussed on increasing the quantity and quality of data, making reciprocal links between ALFRED and other related databases, and providing useful tools to make the data more comprehensible to the end user. ALFRED is accessible from the Kidd Lab home page (http://info.med.yale.edu/genetics/kkidd/) or from ALFRED directly (http://alfred.med.yale.edu/alfred/index.asp).

Molecular markers reveal cryptic sex in the human pathogen Coccidioides immitis.

Coccidioides immitis, cause of a recent epidemic of "Valley fever" in California, is typical of many eukaryotic microbes in that mating and meiosis have yet to be reported, but it is not clear whether sex is truly absent or just cryptic. To find out, we have undertaken a population genetic study using PCR amplification, screening for single-strand conformation polymorphisms, and direct DNA sequencing to find molecular markers with nucleotide-level resolution. Both population genetic and phylogenetic analyses indicate that C. immitis is almost completely recombining. To our knowledge, this study is the first to find molecular evidence for recombination in a fungus for which no sexual stage has yet been described. These results motivate a directed search for mating and meiosis and illustrate the utility of single-strand conformation polymorphism and sequencing with arbitrary primer pairs in molecular population genetics.

Inferring identify from DNA profile evidence.

The controversy over the interpretation of DNA profile evidence in forensic identification can be attributed in part to confusion over the mode(s) of statistical inference appropriate to this setting. Although there has been substantial discussion in the literature of, for example, the role of population genetics issues, few authors have made explicit the inferential framework which underpins their arguments. This lack of clarity has led both to unnecessary debates over ill-posed or inappropriate questions and to the neglect of some issues which can have important consequences. We argue that the mode of statistical inference which seems to underlie the arguments of some authors, based on a hypothesis testing framework, is not appropriate for forensic identification. We propose instead a logically coherent framework in which, for example, the roles both of the population genetics issues and of the nonscientific evidence in a case are incorporated. Our analysis highlights several widely held misconceptions in the DNA profiling debate. For example, the profile frequency is not directly relevant to forensic inference. Further, very small match probabilities may in some settings be consistent with acquittal. Although DNA evidence is typically very strong, our analysis of the coherent approach highlights situations which can arise in practice where alternative methods for assessing DNA evidence may be misleading.

Tsetse thrombin inhibitor: Bloodmeal-induced expression of an anticoagulant in salivary glands and gut tissue of Glossina morsitans morsitans

The tsetse thrombin inhibitor, a potent and specific low molecular mass (3,530 Da) anticoagulant peptide, was purified previously from salivary gland extracts of Glossina morsitans morsitans (Diptera: Glossinidae). A 303-bp coding sequence corresponding to the inhibitor has now been isolated from a tsetse salivary gland cDNA library by using degenerate oligonucleotide probes. The full-length cDNA contains a 26-bp untranslated segment at its 5′ end, followed by a 63-bp sequence corresponding to a putative secretory signal peptide. A 96-bp segment codes for the mature tsetse thrombin inhibitor, whose predicted molecular weight matches that of the purified native protein. Based on its lack of homology to any previously described family of molecules, the tsetse thrombin inhibitor appears to represent a unique class of naturally occurring protease inhibitors. Recombinant tsetse thrombin inhibitor expressed in Escherichia coli and the chemically synthesized peptide are both substantially less active than the purified native protein, suggesting that posttranslational modification(s) may be necessary for optimal inhibitory activity. The tsetse thrombin inhibitor gene, which is present as a single copy in the tsetse genome, is expressed at high levels in salivary glands and midguts of adult tsetse flies, suggesting a possible role for the anticoagulant in both feeding and processing of the bloodmeal.

Nonlinearity in genetic decoding: Homologous DNA replicase genes use alternatives of transcriptional slippage or translational frameshifting

The τ and γ subunits of DNA polymerase III are both encoded by a single gene in Escherichia coli and Thermus thermophilus. γ is two-thirds the size of τ and shares virtually all its amino acid sequence with τ. E. coli and T. thermophilus have evolved very different mechanisms for setting the approximate 1:1 ratio between τ and γ. Both mechanisms put ribosomes into alternate reading frames so that stop codons in the new frame serve to make the smaller γ protein. In E. coli, ≈50% of initiating ribosomes translate the dnaX mRNA conventionally to give τ, but the other 50% shift into the −1 reading frame at a specific site (A AAA AAG) in the mRNA to produce γ. In T. thermophilus ribosomal frameshifting is not required: the dnaX mRNA is a heterogeneous population of molecules with different numbers of A residues arising from transcriptional slippage on a run of nine T residues in the DNA template. Translation of the subpopulation containing nine As (or +/− multiples of three As) yields τ. The rest of the population of mRNAs (containing nine +/− nonmultiples of three As) puts ribosomes into the alternate reading frames to produce the γ protein(s). It is surprising that two rather similar dnaX sequences in E. coli and T. thermophilus lead to very different mechanisms of expression.